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Abstract 

We investigate the spectrum above the kink ground states of the 
spin J ferromagnetic XXZ chain with Ising anisotropy A. Our main 
theorem is that there is a non-vanishing gap above all ground states of 
this model for all values of J. Using a variety of methods, we obtain 
additional information about the magnitude of this gap, about its be- 
havior for large A, about its overall behavior as a function of A and its 
dependence on the ground state, about the scaling of the gap and the 
structure of the low-lying spectrum for large J, and about the existence 
of isolated eigenvalues in the excitation spectrum. By combining in- 
formation obtained by perturbation theory, numerical, and asymptotic 
analysis we arrive at a number of interesting conjectures. The proof 
of the main theorem, as well as some of the numerical results, rely on 
a comparison result with a Solid-on-Solid (SOS) approximation. This 
SOS model itself raises interesting questions in combinatorics, and we 
believe it will prove useful in the study of interfaces in the XXZ model 
in higher dimensions. 
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1. Introduction 

The subject of our paper is the ferromagnetic XXZ model. The XXZ 
model is one of the best studied quantum spin systems, benefiting from 
both algebraic and analytic techniques. But, while the mathematical 
techniques that have been applied to the XXZ model are impressive (c.f. 
in many of the most basic physical questions remain open ||22|| . 



We address an unresolved issue of the one-dimensional XXZ model, 
which is the following. It is known, by rigorous methods |ToT|, that 



there is a spectral gap above the infinite volume ground state for spin 
1/2, but the proof relies on an algebraic tool which is not present for 
higher spins. How does one prove the existence of a spectral gap in the 
more general setting? We answer the question in the present paper. 
Our methods are somewhat more general than those of ]16[ since we 
do not rely on the quantum group symmetry. On the other hand, it 
is essential for our proof that we know the spectral gap exists for the 
spin 1/2 XXZ model. Still, we believe our techniques may be applied 
to other spin models, as well as shedding light on this corner of the 
general knowledge of the XXZ model. 

The XXZ spin chain is a generalization of the Heisenberg model 
where one allows anisotropic spin couplings. The Hamiltonian for the 
spin J model is 

(i-i) #i J) = - £ o&sJ + s « s l + AS « S D > 

where S^ 2 ' 3 are the spin J matrices acting on the site a, tensored with 
the identity operator acting on the other sites, (a, (3) denotes a pair 
of nearest neighbors. The local Hilbert space is 7i a = C 2J+1 , and 
is a Hermitian operator on H\ = (SWa^q- ^ or now we think of A 
as a finite subset of Z, though we are also interested in the case that 
A = N or Z. The main parameter of the model is the anisotropy 
A e R. By choosing A = ±1 we can obtain the isotropic ferromagnet 
or antiferromagnet. Alternatively, by taking A — > ±oo we recover 
the spin J Ising ferromagnet and antiferromagnet. In this paper we 
restrict A > 1, which corresponds to a ferromagnet with the strongest 
coupling along the S^-axis. We note that H\(A) and — H^(— A) are 
unitarily equivalent. It is useful to introduce two other forms of A: 
q = A — VA 2 — 1 and r\ = — logg. Observe that A = |(g + g" 1 ) = 
cosh(?7). The parameter q is the one which labels the quantum group 
SU g (2) when J = 1/2. As A increases from 1 to oo, q decreases from 1 
to 0, and r\ increases from to oo. 
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For A > 1, it is widely known that there are two infinite volume 
ground states which correspond to all spins up, |+J) and all spins 
down, |— J). It is considerably less well known that there are in fact 
many more infinite volume ground states. These extra ground states 
come in two families: the kink states and the antikink states. The 
kink states are an infinite family of ground states all with the same 
GNS space, which break discrete translation symmetry as well as the 
continuous U(l) symmetry associated to the XXZ model. They have 
the property of being asymptotically all down spins at — oo and all 
up spins at +00. They clearly also break left-right symmetry: their 
reflected counterparts are the antikink states. Our main results concern 
the spectral gap above these kink ground states. The kink states are 
physically interesting for several reasons: They exhibit domain walls, 
an important feature of real ferromagnets (see |25j for an application of 
the XXZ model to spin droplets in one dimension); For J > 1 the kink 
ground states of the XXZ model are more stable than the Ising ground 
states, which is a new result and subject of the present paper. Also, the 
XXZ spin chain plays an important role in explaining the phenomenon 
of negative resistance jumps and hysteresis in recent magnetoresistance 



experiments Jl8|, j29j. 



Although our main subject is the XXZ model for J > 1/2, let us 
briefly recall some important facts about the spin 1/2 model. What is 
probably most well known is that the spin 1/2 model is Bethe ansatz 
solvable. This is not applicable to our case however, and we will not 
use the Bethe ansatz in any way. A second interesting feature of the 
spin 1/2 XXZ model is that it possesses a quantum group symmetry. 
Specifically, in PBf it was shown that adding a boundary field 



(1.2) B = JVA 2 - 1(5 X 3 - Si) 

makes the Hamiltonian commute with SU g (2) on 7i\. (Actually there 
are two representation of SU g (2) corresponding to the two opposite 
linear orderings of its tensor factors; H\ + commutes with one 
representation and H\ — B\ commutes with the other. We will only 
consider H\ + B\.) For < q < 1 the representation theory of SU g (2) 
is equivalent to that of SU(2) (c.f. ]15|), and it plays the same role in 
the analysis of the XXZ model that SU(2) plays in the analysis of the 
isotropic model. For example, the ground state space corresponds to 



the highest- dimensional irreducible representation of SU 9 (2). In [16 
Koma and Nachtergaele used the quantum group symmetry to calculate 
the spectral gap for the XXZ model by proving that the lowest exci- 
tations of the XXZ model form a next-highest dimensional irreducible 
representation of SU 9 (2). There are still open conjectures relating to 



4 



the representations of SU 9 (2) and the XXZ model such as: Prove the 
lowest (highest) energy of H\ + B\ restricted to the spin s representa- 
tions in 7Ya is lower than the lowest (highest) energy of the spin s — 1 
representations. This is almost certainly true, and would generalize the 
Lieb and Mattis result [2D[, but remains open. 

We now turn our attention to the spin > 1/2 models. The first im- 
portant fact is that the ground states have been explicitly calculated 
for all finite volumes, in all dimensions, all choices of spin, and even 
allowing different values of anisotropy along bonds in the different co- 
ordinate directions UJ. In |T2[ the ground states were independently 
dicovered for J = 1/2 and one dimension, and these ground states were 
generalized to infinite volume ground states. They have the property 
of being frustration free, which means that they not only minimize 
the expectation of the infinite volume Hamiltonian, they minimize ev- 
ery nearest neighbor interaction, too. Gottstein and Werner found all 
the frustration free ground states, and conjectured that there were no 



other ground states. In [21|, their conjecture was proved correct, and 



the analogous statement for J > 1/2 was proved in [17]. Thus, one has 
a complete list of ground states for the one dimensional XXZ model 
with any choice of J. Unfortunately, the results on infinite volume 
ground states are valid only in one dimension. Finding the complete 
set of ground states in dimensions two and higher is an important open 
problem. 

In this paper we prove that there is a nonvanishing spectral gap 
above the infinite volume ground states for every J, thus extending 
the results of |16|]. We mention that the existence of a spectral gap is 



generally believed to follow from the fact that the quantum interface 
of the kink ground states is exponentially localized. Our results verify 
the conventional wisdom, and our proof does rely on the exponentially 
localized interface. However there are other important elements to our 
proof: most notably, a rigorous comparison of the spin J chain with 
a spin 1/2 ladder with 2 J legs. The spin J XXZ chain is a quantum 
many body Hamiltonian. The Hilbert space for the L-site spin chain 
is (2J + 1) L , its dimension grows exponentially with L. The dimension 
of the spin ladder is even larger, at 2 2XL . But an important bound, 
Lemma |4~l"l , allows us to restrict attention to an (L + 1) 2J dimensional 
subspace. This allows the proof of the existence of the spectral gap, 
and also allows more efficient numerical methods for studying the XXZ 
model. The reduced system resembles a quantum solid on solid model 
for the spin ladder. We view the present problem as a warm up for 
the QSOS method, which we believe will play an important role in 
proving stability of the 111 interface for the XXZ model. We also 
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expect the spin ladder technique will be useful in proving the existence 
of a spectral gap for other spin chains where a gap is known in J = 1/2 
but not for J > 1/2. In this paper, in addition to giving a rigorous proof 
of the existence of a gap, we present a new type of numerical method 
for studying the XXZ model. We also present an asymptotic model for 
the low lying spectrum of the XXZ model as J — > oo, in terms of a 
free Bose gas. The asymptotics explain new qualitative features of the 
XXZ model for J > 1/2, and is in excellent agreement with numerical 
data for J sufficiently high. 

The remainder of the paper is organized as follows. In Section 2 we 
present our main theorem, as well as a number of conjectures which are 
supported by numerical evidence and asymptotic analysis. In Section 
3 we introduce some background material which is useful for our proof. 
In Section 4 we derive the spin chain / spin ladder reduction. In Section 
5 we finish the proof of the main theorem. In Section 6 we combine the 
lower bounds for the spectral gap with numerical methods to obtain 
data for the spectral gap. In Section 7 we derive a boson model for 
the XXZ spin system which explains the asymptotic behavior of the 
gap as J — > oo. This boson model is similar to [13|, ||, [K|, but without 
the need for a large external field (other than the boundary field which 
vanishes in the thermodynamic limit). 

2. Main Result and Conjectures 

(The notation [a, b] will always refer to the discrete interval {a, a + 
1, . . . , b}. It is not necessary that a and b are integers as long as the 
difference b — a is.) 

The Hamiltonian we will use is the following spin J XXZ Hamilton- 
ian: 

{a,«+l}cA 

( 2 - 3 ) h\a, a + 1) = (j 2 - S 3 a S 3 a+l - A^SX+i + S 2 a S 2 a+1 ) 

+ iVT^^(sl - s 3 a+1 )) . 

In comparison to the Hamiltonian fll.l] ), we have just added the bound- 
ary fields ( |1.2|) , scaled by A" 1 and added a constant. One can easily 
check that the interaction h\a,a + 1) is nonnegative. For finite vol- 
ume A, it is an easy but important observation that the Hamiltonian 
commutes with S\ = YlaeA^a- We use this symmetry to block diag- 
onalize H^. In particular, we let 7i(A, J) be the spin J Hilbert space, 
and we define TC(A, J, M) to be the eigenspace of S\ with eigenvalue 
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M £ {— J|A|, . . . , J|A|}. We call these subspaces "sectors". They are 
invariant subspaces for H^. 

For finite volumes A, the ground states of Hj^ may be expressed in 
closed form, as was pointed out in [jlj . We will give a formula for these 
ground states in the next section. For now we merely mention the fact 
that for each sector there is a unique ground state \I/o(A, J, M), and it 
has the property that its energy is zero. We define the ground state 
space Q(A,J,M) to be the one-dimensional span of ^o(A, J, M), then 
the spectral gap is given by 

(2.4) 7(A,J,M)= mf 

</>eW(A,j,M) 

«/>_Lg(A,J,M) 

where (• • = (^|- • 

One passes to the thermodynamic limit, by considering the infinite 
volume Hamiltonian as the generator of the Heisenberg dynamics on 
the algebra of quasilocal observables Aq. The definition of a ground 
state is a state on Aq such that for any local observable X £ Aa, 
|A| < oo, uj satisfies 

(2.5) uj(X*5(X)) > 

where 5{X) = lim\yz[H^, X]. As in the case of finite volumes, there 
is a collection of ground states whose GNS representation is explicit. 



These ground states were discovered in [12], and they were proven to 



be the complete list in [|17]]. The infinite volume ground states are 
the following: a translation invariant up spin state determined by the 
equation uj^(S^) = +J for all a; a translation invariant down spin state 
ujt; an infinite number of kink states which we label u;^; and an infinite 
number of antikink states, uj\\. The kink states have the property that, 
if T is the translation to the left one unit (so T~ 1 S^T = S^ +1 ), then 
for any quasilocal observable X 

lim uH(T- n XT n ) = uj\X), 

(2.6) ™ .. 

lim cj|j(T n XT~ n ) =cu l (X). 

n— >oo 

The label M is any integer and is determined as follows. All the kink 
states are also local perturbations of one another, which we will see in 
the next section when we write the explicit GNS representation. So 
there is only one GNS Hilbert space for all the kink states, and only 
one GNS Hilbert space for all the antikink states. 

For any ground state u, the infinite volume Hamiltonian can be 
represented as the generator of the Heisenberg dynamics for the algebra 
of observables on Hqns, the GNS Hilbert space of uj. This means that 
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there is a densely defined, self adjoint operator T/gns with the property 
that for any X € Ao, 

#GNS7TP0^GNS = 7t(^P0)^GNS , 

where Ogns is the representation of the ground state as a vector, n is the 
representation of the quasilocal observable algebra on the observable 
algebra of 7Ygns, and S(X) = \im^y"z[H^,X] is the derivation defining 
the Heisenberg dynamics. The bottom of the spectrum of -£/gns is 0. 
The spectral gap above u> is defined to be the gap (if one exists) above 
in the spectrum of H GNS . Note that there is one spectral gap for 
each of the four classes of ground states: all up, all down, kinks, and 
antikinks. 



We can now state the main result of [IB 



Theorem 2.1. (|T(J For the SU q (2) invariant spin-1/2 ferromagnetic 
XXZ chain with the length L > 2 and A > 1, the spectral gap is 

7([1, L], 1/2, M) = 1 — A" 1 cos(tt/L) , 

in any sector H([l, L], 1/2, M), -L/2 < M < L/2. Above any of the 
infinite-volume ground states (all up, all down, kink or antikink) the 
spectral gap is 

7 = 1 - A" 1 . 

■ 

The previous theorem is the starting point of our own analysis. Our 
main result is an analogous theorem, extending the existence of the 
spectral gap to all J instead of just J = 1/2. 

Theorem 2.2. For any J G ~N ; and any A > 1, the gap above the 
translation invariant ground states is 

lup = Idown = 2J(1 - A _1 ) . 

The gap above the kink states satisfies the bounds 

^ Ikink — 7«p • 

Specifically the spectral gap above the kink ground state is nonvanishing. 

Remark In the theorem, the formula for 7 up is well-known, The in- 
equality 7kink < 7u P is also well-known and easy to deduces. We include 
proofs of these facts for the convenience of the reader. Our main result, 
which is new, is that 7km is strictly positive. 

We can say something more specific about the low excitation spec- 
trum by considering an extra symmetry of the Hamiltonian. Since Hj^ 
commutes with S\ for each finite volume A, we would like to define 
an infinite volume analogue of S\. For the GNS space above the kink 
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ground states the correct definition is the following renormalized ver- 
sion 

which is a densely defined self adjoint operator on the TCqns- The 
ground state space of -£/gns is spanned by the orthogonal family of 
vectors {^>q(1j, J, M) : M G Z} which are determined up to scalar 
multiplication by the properties that 

# GNS *o(Z, J, M) = , S$V (Z, J, M) = MV (Z, J, M) . 

We define a version of the spectral gap for Hamiltonian restricted to 
the sectors of S 3 in the following way. Let 7(Z, J, M) be the largest 
number such that for any local observable X G Aa commuting with 
S\ we have 

(*o(Z, J, M)\ir(X)*H 3 GNS ir(X)y (Z, J, M)) 

> 7 (Z, J, M)<* (Z, J, M)|7r(X)*^ NS 7r(X)* (Z, J, M)> . 

An arbitrary local observable does not commute with S\. However, one 
may define Xm for — J|A| < M < J|A| so that each Xm commutes with 
Si and vr(X)^o(Z, J, M) = n{X M >)Vo{Z, J, M + M'). This is just 
due to the fact that the GNS representation (7Ygns, tTj ^o(Z, J, M)) is 
cyclic for any choice of M. From this we see that 

7kink = mf^7(Z, J,M) . 

Let T be translation to the left, as before. We have T _1 S 3 T = 
2 J + S 3 , which implies 

T^ (Z, J, M) = * (Z, J, M + 2J) , 

since T clearly commutes with the Hamiltonian. Hence, 

7 (Z,J,M) = 7 (Z,J,M + 2J). 

Another symmetry of the Hamiltonian is obtained by taking a left-right 
reflection of the lattice about the origin, and simultaneously flipping 
the spin at every site. This is a unitary transformation of 7Ygns to 
itself. Calling this symmetry TZ we have 1ZS 3 1Z = — S 3 . So 

7 (Z,J,M) = 7 (Z,J,-M). 

To prove that the gap above the infinite volume kink states is nonzero, 
it suffices to check that 



7(Z, J,M) > 
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Figure 1. The spectrum of the XXZ spin chain for 
various finite spin chains. The horizontal axis is the 
total magnetization fif ot of a sector: the lines above 
each number show the eigenvalues in that sector. The 
spin systems are, from left to right and top to bottom, 
(J, L, A) = (1/2, 10, 2); (1, 7, 4); (3/2, 6, 4); (2, 5, 8). 



for M = 0, 1, . . . , [J] , where [J] is the least integer greater than J. This 
fact does not actually simplify the proof, but the symmetries above are 
an important part of our proof. 

In Figure [l] we show the spectrum for some small spin chains, as 
calculated by Lanczos iteration. Even though the lengths are finite, 
one sees that the gap is an even function of M, and is nearly periodic 
of period 2J. As one takes L — > oo, the gap for any finite region of M 
values is periodic and even. 
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Figure 2. Plot of 1, 1-5 (left) and 0, 7 (right) ver- 
sus A^ 1 , for (J,n, L) equal to: First column (7/2,1,4), 
(2,0,6), (4,1,4), (5/2,0,5); Second column (3,0,4), 
(7/2,0,4), (4,0,4), (9/2,0,3) 

Our main theorem proves existence of a spectral gap, but it is obvi- 
ously just as interesting to know what the gap is. Unfortunately, the 
most information we can gain from our proof is that the spectral gap 
can be well approximated by calculating the gap in a finite volume, L, 
with an error which decreases like q L . This still leaves the problem of 
calculating the gap in a finite volume L oc 1/rj. We do not have any 
rigorous bounds for the spectral gap valid for all g £ (0, 1). However we 
have studied the problem in three ways: numerically, by perturbation 
series, and asymptotically; and we propose the following conjectures 
based on our findings. 



2.1. Numerical results. We performed two types of numerical meth- 
ods. The first, and more efficient method is based on the spin ladder 



reduction from our proof. In Theorem |4.1| below, we obtain a rigorous 
lower bound for 7QI, L], J, M) as 



7 ([1, L],J, M) > 2J(1 - A^)(l — S([l, L], J, M)) 




Figure 3 . Some plots of the spectral gap using Lanczos 
iteration. For all three curves, n = mod 2J. 

where 1 — 5([1,L],J,M) is the spectral gap for a reduced model re- 
sembling a quantum solid-on-solid model for the spin ladder. What is 
important about the bound is that the reduced model has dimension 
(L + 1) 2J as opposed to the original system with dimension (2J + 1) L . 
We then numerically diagonalized the reduced system to find the spec- 
tral gap for some values of J and L. The second numerical method was 
simply to numerically diagonalize the original Hamiltonian for some 
small values of L and J. We did this primarily to check the qualitative 
results of the lower bound. In Figure |2] we show the results of the lower 
bound calculation, and in Figure |3| the result of the Lanczos iteration. 
What emerges qualitatively is that for J > 1 there is a local maximum 
for the spectral gap with 1 < A < oo. This is other than expected 
based on the spin 1/2 results. Based on our numerical evidence we 
make the following conjecture. 

Conjecture 2.3. We have defined 7(Z, J,M) above for fixed A. Let 
us rewrite this as 7(Z, J, M, A~ x ) to take account of the anisotropy 
< A -1 < 1. We conjecture that 

min7(Z, J, M, A" 1 ) = 7 (Z, J, 0, A" 1 ) , 

and that for each J there exists a A7 1 such that 

7(Z,J,0,A- 1 )< 7 (Z,J,0,A7 1 ) 
whenever < A -1 < 1, equality holding only if A -1 = AT . For 



J — 1/2, the conjecture is a known fact following from Proposition \2. 1\ , 
and one sees AT^ = 0. We conjecture that Ar 1 = as well, but that, 
for J > 3/2, < A7 1 < 1. 
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2.2. Ising Perturbation. If we write the Hamiltonian as an expan- 
sion in A -1 , we can make a perturbation expansion off of the Ising 
model. The Ising model spectrum is well known, but some interesting 
facts arise. One fact which is useful is that the Hamiltonian obtained by 
changing A -1 to — A -1 is unitarily equivalent to the original: just ro- 
tate every other site by ir about S 3 . This means that A" 1 = is either 
a local minimum or a local maximum of 7(Z, J, M, A -1 ). For J = 1/2, 
the first excitation above a kink ground state in the Ising limit is infin- 
itely degenerate. This is why the slope of 7(Z, 1/2, M, A" 1 ) = 1 — A -1 
is not zero at A" 1 = 0. Similarly for J = 1 and M odd. However, for 
all other choices of J and M, the first excitations are at most finitely 
degenerate, and for J > 1 and M = the first excitations are nondegen- 
erate. This means that the first derivative of 7(Z, J, M, A" 1 ) vanishes 
for all other values of J and M. For J = 1, M = the second derivative 
is negative, while for J > 1 and M = the second derivative is always 
positive, indicating that the Ising limit does not maximize the spectral 
gap, but minimizes it locally. This argument, which will be expanded 



in Section 6, illuminates part of Conjecture |2.3| . The nondegeneracy 
of the first excitations suggests a second gap. Based in part on this 
evidence, we make the following conjecture: 

Conjecture 2.4. 1) For l<A<oo ; J > 3/2 and any MeZ, the 
lowest excited state is an isolated eigenvalue, i.e. there is a nonvanish- 
ing gap to the rest of the spectrum. 

2) For J = 1 and M any odd integer, the lowest excited state is the 
bottom of a branch of continuous spectrum. For J = 1 and M even, 
the lowest excited state is again an isolated eigenvalue. 

2.3. Asymptotics. From the numerics it became clear that for M = 
and A fixed, the gap 7(Z, J, 0, A -1 ) scales like J. Also, a careful 
analysis of the exact formula for the ground states of H, which will 
be presented in the next section, shows that for large J, the wave 
vector has Gaussian fluctuation on the order of J 1 / 2 . These two facts 
together suggest a scaling analysis of the bottom of the spectrum of 
Hj^ in the limit J — > oo. Consistent with the Gaussian form for the 
ground states, our asymptotic analysis leads to a free Boson gas model 
for the bottom of the spectrum, at least to first order in J~ 1 / 2 . We 
derive a boson model for the XXZ spin system analogous to (See 
also 0, 10] for a better introduction to spin waves. Unfortunately our 



treatment is not as well developed.) This is not the same analysis 



as was done in ||19|| , nor in any other coherent states approach. We 
analyze the eigenstates whose energy scales like J, whereas coherent 
states give rigorous bounds on the bulk spectrum which scales like 
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J 2 . The higher energy states are much greater in number, so typically 
they control the thermodynamic behavior. However, the low energy 
states may be more important for dynamical properties (cf ||), since 
that part of the spectrum is separated by spectral gaps. Also note, the 
asymptotic model is that of a free Bose gas, with a nontrivial dispersion 
relation for the energy of each oscillator. In fact, the energy is such that 
excitations which give least energy are exponentially localized about 
the interface. We believe that this makes it plausible to prove the Boson 
gas estimate is correct with small errors for large but fixed J < oo, 
because even as one takes L — > oo, the low excitations are "essentially 
finite" and we can more or less prove the large J asymptotics for the 
finite system. Of course J would have to depend on the number of 
excitations that you wanted to estimate by the Bose gas picture; as 
the number of excitations goes to oo so must J. Our boson model has 
a quadratic coupling, but a Bogoliubov transformation diagonalizes 
it. The spectrum of the coupling matrix, then gives the value for the 
limiting curve of the spectral gap, and other information about the low 
spectrum. Based on our analysis we make the following conjecture: 

Conjecture 2.5. Q Th ere is a function 7oo : (1, 00) xl->R with the 
property that 

lim r^(Z, J, tiJ, A) = 7oo (/i, A) . 

J— >oo 

This function satisfies 

7oo(/i + 2, A) = 7oo(/i, A) and 7<x> (-/■*, A) = 7 00 (/i, A) . 

Moreover 7oo (/i, A) is equal to the spectral gap of a bi-infinite Jacobi 
operator A (really A(A,fi)) defined on Z 2 (Z) 

Ae n = [2e n - sech(?7)e n „i - sech(77)e„ + i] 

(2.7) 4sinh 2 (r/) 

cosh(2?7(n — r)) + cosh(2?7) n 

Above, r = r(/i, A) is a phase defined implicitly by the equation 

n 

\i = lim \ tanh(?7(n — r)) . 

-n+1 

Remark 1. The implicit formula for 7oo (/i, A) has a specific conse- 
quence that 

lim Af 1 = 0.49585399 ± 1(T 8 , 

1 Note Added in Proof: Since the submission of this paper to the arXiv, 
Caputo and Martinelli have obrtained further results which prove part of Conjecture 
2.5. In M they obtain a lower bound for the gap of the form J x constant. 
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J -> °° spectral gap, y , as a function of anisotropy and filling factor 



1 




Figure 4. Surface and contour plot for the function 

7oo versus r and A -1 

which is obtained by numerical diagonalization of large Jacobi matrices. 

Remark 2. Equation ( j2.7| ) has a simple physical interpretation. The 
bracketed term on the RHS is the usual matrix for a one-magnon spin 
wave, which by itself would give the energy 2(1— sech?]) for the spectral 
gap, just as in the translation invariant ground states. The second 
term on the RHS of (|2.7| ) represents the attractive potential due to the 
domain wall with center r. 

In Figure [| we show the function 700(7", A -1 ) as obtained by numer- 
ically diagonalizing the Jacobi operator for 50 sites. Note that the 
number of sites reflects the extreme simplicity of the Bose model over 
the true XXZ model; we could not numerically diagonalize a spin chain 
of 50 sites even for spin 1/2. 

3. Ground states of the XXZ model 

In this section we give formulas for the ground states of the finite 
volume XXZ model as in []]]. We show how the ground states of the 
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spin J chain can be derived by looking at the ground states of a spin 
1/2 ladder with 2 J legs, which is essential for our proof. We also recall 
the explicit GNS representation of the infinite volume ground states in 
a Guichardet Hilbert space (also called incomplete tensor product) [12|, 



1 
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17j|. This concrete representation is convenient, especially for proving 
that certain sequences of finite volume ground states have unique limits. 

We begin by rewriting ( |2.3| ) for spin 1/2. It is easy to check that for 
the two site interaction h}/ 2 (l, 2), the three vectors 

l + 2 >+ 2 )> I 2 > 2 ), I 2 >+ 2 ) + ?l + 2? 2 ) 

are ground states, while |+l/2, — 1/2) — q\— 1/2, +1/2) is a state with 
energy one. In other words, we may view /^(M) as I — U(t(1, 2)), 
where t(1,2) G ©2 is the transposition, and U is the (non-unitary) 
action of 62 defined by 

f/(r(l,2))0(mi,m 2 ) = 0(m 2 ,mi), 
where the 0(m!,m 2 ) are the (non-normalized) basis vectors 

-i>, «-i,i) = 8 v.|4i>. 

In other words, the ground states of the two site Hamiltonian are the 
symmetric tensors with respect to the nonunitary action U. We can 
generalize this result to linear chains of any length, and to many other 
domains as well. Specifically, what we need to properly define the 
XXZ Hamiltonian with boundary fields is a collection of sites A and 
a collection of oriented bonds among those sites B, in other words a 
digraph. Then we define 

H{,b= E h * 

(a,/3)eS 

(We write when B is obvious.) We define a height function to be 
any function I : A — > Z such that Z(/3) — Z(a) = 1 for all (a, (3) G B. 
The condition to have such a height function is that for any closed 
loop, where a loop is defined as a sequence a±, at2, ■ ■ ■ , a n = ol\ G A 
such that for each i either (ati, ccj+i) G £> or (a,j + i,aj) G £>, there are 
equal numbers of bonds with positive orientation (aj,a i+ i) e B as 
with negative orientation (aj+i, ai) G B. The following lemma is an 
interpretation of a result in M. 
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Lemma 3.1. If (A, 23) is a connected digraph such that a height func- 

1/2 

tion I exists, then there is a unique ground state of H^ B in each sector 
W(A, 1/2, M) for M G [— |A|/2, |A|/2]. 

Proof: The proof is like the analogous statement (without the re- 
quirement of a height function) for the isotropic model. Suppose that 
/ exists. Define a (non-normalized) basis of vectors 

<f>({™«}) = q-^ Amalia) \{m a }}. 

Then define an action of ©a on H.(A, 1/2, M) by 

U(n)(f>({m a }) = 0({m ff -i( a )}) . 

As we have already seen, for any (a, (3) G 23, h l / 2 {a,(3) = I— 27(r(a, /?)), 
where r(a,P) is the transposition. Thus, the unique ground state vec- 
tors are those vectors which are symmetric under the action of U(&\), 

tt (A,l/2,M)= £ 0({m Q }), 

{m a }e[-l/2,l/2] A 
J2 a m a =M 

one for each sector. ■ 

We note that one can trivially prove the converse of this lemma, that 
if there exist ground state vectors in any sector other than M = ±|A|/2, 
then there is a height function (as long as < q < 1). We now 
mention a second lemma (which is also implicit in [E1J) which gives the 
construction of the spin J ground states by using spin ladders. 

Lemma 3.2. Suppose that (A, 23) satisfies the hypotheses of the last 
lemma. Then for any spin J G |N ; there is a unique ground state of 
Hjlg in each sector T~C(A, J,M). Moreover, this ground state is associ- 
ated to the ground state of the spin 1/2 spin ladder 

A = {(a J) : a G A,j G [1,2 J\}, 

B = {((a, j), ((3, k)) : (a, (3) G 23, (j, k) G [1, 2 J] 2 } . 

Define Q a : H({a} x [1,2.7], 1/2) — > H({a}, J) to be the projection 
onto the unique highest spin representation in the decomposition of 
7i({a} x [1,2 J], 1/2) into irreducibles, so that Q* a Q a is the projection 
onto symmetric tensors. Then Q\ = n^eA Qa 9^ ves an isomorphism 
of ground states 

*o(A, J, M) = g A ^ (A, 1/2, M) , * (A, 1/2, M) = Q* A V (A, J, M) . 



17 



Proof: Most of the proof is self evident, the point being just to 
introduce the notation necessary for later work. We note that we can 
define a height function I by l(a,j) = 1(a). Then the ground states of 

1/2 

H~g are in the range of Q\Qk because for any permutation it which 

preserves the rungs of the ladder, i.e. vr({a} x [1,2J]) = {a} x [1,2J] 
for all a, U(ir) coincides with the normal action. So 

(3.1) Q* A Q A V (A, 1/2, M) = * (A, 1/2, M) • 

On the other hand Q A H^ 2 Q* A equals H{. So Q A *o(A, 1/2, M) is a 
ground state of H A , in fact all the ground states are obtained like this by 
an easy Perron- Frobenius argument. Rewriting ( |3.1| ) with ^o(A, J, M) 
in place of Q\^q(A, 1/2, M) finishes the proof. | 



In the particular case of A C Z, one can take 1(a) = a. Let us for 
future notational ease define Ai(A, J, M) to be the set of all {m a } £ 
[— J, J] A with the property that X] a m a = M. Then the formula one 
derives for the spin J ground states is 

* (A,J,M) = E II (j+ J m ) V \— |K}>. 

We now turn our attention to the infinite volume ground states. We 
will merely define the ground states, the proof of completeness was 
done in fll7f . Given a countably infinite set of sites Aoo, and a finite 
dimensional Hilbert space 7i a and unit vector Q a at each site, one can 
define the Guichardet Hilbert space 

(oo n oo 

n=l j=l j=n+l 
n— 1 oo 

n(0w« 3 . ®(g)Cfi a ,) 

j=l j= n 

where a±, ■ ■ ■ is any enumeration of A ro , and cl means the usual L 2 
closure. If A is any finite subset of A^ then we can define the finite 
dimensional Hilbert space Ha = ®«eA^a' as usua l ; an d an obvious 
inclusion 

ia,a x ■■ n A ^ (g) (n a ,n a ). 

This is the proper framework to discuss the one dimensional infinite- 
volume ground states of the XXZ model, because of the following 



18 

Definition 3.3. For any J £ |N we make the following definitions. 
Consider = N. For each a £ Zet 7Y Q = C 2J+1 and Q a = \J). 
Denote the Guichardet Hilbert space so obtained by 

w(n, j,«p) = (g)(w a ,n a ). 

Now let Aqo = Z instead. For each a £ Aoo let Ti a = C 2,/+1 as before, 
but let Q a = | J) for a > 1 and f2 Q = | — J) /or a < 0. Define 

H{Z,J,kmk) = 0(ft Q ,n a ). 

Then 

Lemma 3.4. (a,) // a sequence of integers with L k — > oo, and 

A" £ N, £/ien £/ie normalized sequence 

i[i,L k ]^o([i-,L k ],J,JL k -N) 
\\M[hL k ],J,JL k -N)\\ 

converges in norm. 

(b) If M k = L k J- 2Jr k + N where N £ [0, J - 1], r k £ N and both r k 
and L k — r k tend to oo, then the normalized sequence 

T- r H [hW V ([l,L k },J,M k ) 

\\* ([l,L k ],J,M k )\\ 

converges in norm, where T is the translation one unit to the left. 

Proof: We prove this for spin 1/2. Then the analogue follows for 

1/2 

ground states of H~ , and by the last lemma this proves it for arbitrary 
spin. 

For (a), note that defining 

L], 1/2, N) — 9 ai+ - +ow ^ " " " S~ N n N , 

l<ai<--<a N <L 

where is the all up spin vector, we have 

» [1|L] , z tto([l,L],l/2,±L-iV) ^([1,L],1/2,AQ 
\\M[1,L],1/2,±L-N)\\ \\%([1,L],1/2,N)W 

Then (a) is proved if we prove that the sequence ^([1, L], 1/2, N) 
converges. But this follows by the Monotone Convergence Theorem, 
thinking of the coefficient of S ai ■ ■ ■ S aN Q^ as a function ■ ■ ■ ,c*n)- 

We still need to check that the limit is finite, i.e. that 

q 2(a 1+ -+a N ) < ^ _ 

l<ai<--<ajv 



1!) 



We can evaluate the series explicitly; it is q N ( N+l ^ / Y{J =1 (1 — g 2j ). 

For (b), we do a similar thing. We define Q% to be the vector 
<S) a€ z Since J = V 2 , now M k = \L k - r k where r k , L k -r k ^ oo. 
Then 

T- r H [hLk] ^ ([l,L k ],l/2,M k ) V' Q ([l-r k ,L k -r k ], 1/2,0) 
||tt ([l, L k ], 1/2, M fc )|| " - r k , L k - r k ], 1/2, 0)|| 

where 

oo n 

%([-a : 6], 1/2,0) = £ £ n^"° fc ^^ n *- 

n=0 -a<ai< --<a„<0</3i< --</3„<6 fe=l 

The lemma will follow by the Dominated Convergence Theorem if we 
prove that ^q(Z, 1/2,0) is summable. This is equivalent to 
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Z = ^ ^2 q 2 lPi+"+Pn-(<*i+~+a n )] < ^ 

n=0 Qi<--<a n <0</3i<-< / 8 n 

We can also evaluate this explicitly 

^n(n+l) qn(n— 1) 



7 - V * 

^n- =1 (i 



1 



by Heine's theorem (c.f. f lj). In particular it is finite. 



Let 



i-=o ||*„([-L+1,L],J,M)|| ■ 
The limit exists by the lemma. We claim that this vector gives an 
infinite volume ground state. The representation of quasi-local ob- 
servables on 7i(Z, J, kink) is clear: the local observables A\ are op- 
erators of 7Ya which includes by into TC(Z, J, kink). Take the 
weak-* completion and we are done. The state is then ujf^iX) = 
(*o(Z, J,M)\XV (Z, J,M)). All one needs to check is that for any 
local observable X, 

lim (V (Z,J,M)\X*[H A ,X]V (Z,J,M)) > 

A^oo 

Suppose X e ^4a- Define Ai to be the union of A with all nearest neigh- 
bors. Then for A' D Ai, [H A /,X] = [H Al ,X]. Using the frustration free 
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property of \&o([ — L, L], J, M), we see that as soon as A 2 C [— L, L], we 
have 

(* ([-L+l,L],J,M)\X*[H Aa ,X\V Q ([-L + l,L],J,M))>0. 

Taking the limit, we see that ^j™/ is an infinite volume ground state. 
We define the r y{ J, M) to be the gap above zero in the spectrum of the 
Hamiltonian acting on the GNS space of the ground state ojf^. 

There are three other classes of ground states. The antikink ground 
states are the states ujf^ = uj^m ° 3~ where T is uniquely determined 
by the formula T{S^ = for all fceZ. There are also the all up spin 
states and all down spin states, which are well known and characterized 
by Lu^ p ' down (Sl) = ±J for all a. We mention that the GNS space for 
them is also a Guichardet Hilbert space where Q™p> down = \±J) for all 
a, and f2^ P ' ° wn is the vector representing the ground state. That these 



are all the ground states is proved by Koma and Nachtergaele [17j 



4. Spin ladder reduction 

We now elaborate on the spin ladder construction introduced in 
Lemma |3.2| . Let us define H-^ = H-^ g introduced in the lemma. Then, 
defining (?(A, J, M) to be the one dimensional ground state space of 
H{ in the sector H(A,J,M), and G(A,M) to be the one dimensional 
ground state space of 7~C(A, M), the lemma tells us that 

Q A Q(A, M) = g(A, M) , Q* A G(A, M) = £(A, M) = Q(A, M) . 

The spectral gap in the sector TC(A, J, M) is defined 

7 (A,J,M)= inf <^M. 

il)±g(A,J,M) 

In view of the lemma, defining Pa = QXQa, we can rewrite the spectral 
gap 



7 (A, J, M) = _inf 



(P A ^\H~ A P^) 



V.eW(A,M)\kerP A (PAlp\ Pa^P) 
jI)±G(A,M) 

We now introduce a second Hamiltonian on Ti.(A, M) which is 

A A,B ' 

where 

&={((a,j),((3,j)):(a,(3)eB, j'e[l,2J]}. 
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This is clearly equivalent to 2 J disjoint copies of H 1 / 2 . So, if A = 

[1, L], Theorem guarantees that the spectral gap of H A is equal to 
1 — A~ 1 cos(7r/L) > 1 — A -1 . Let us introduce some notation: Let 

7Y (A, M) be the ground state space of H A in the sector H(A, M), and 

let 7i cxc (A,M) be its orthogonal complement in the sector H(A,M). 
Then 

(4.2) inf. <*!M>i-A->. 

Ve«cxc(A,M) wm 

Defining P a = Q* a Q a to be symmetrization in the rung {a} x [1,2 J], 
we observe that 

1 2J 

Pa& (a,j)Pa = ^"J ^ ^ Pa& (a,k)Pa ■ 
k=l 

From this it follows that 

2J 



P^h^daJ), (P,j))P a Pp = P a P^ 2 ((a,k), (P,l))P a Pp 

^ ' k,l=l 

and finally that 

P A H~ A P A = ^PaH~ a Pa ■ 
Then we can rewrite the spectral gap formula once more 



(4.3) 7 (A, J, M) = 2J _inf 



(P A ^\H~ A P A ^) 



i/>eH(A,M)\ker P A (Pa^|Pa^) 
ip±g(A,M) 

This is useful because we have a spectral gap for H A . Also note that 
it is now trivial that Q(A, M) C H (A, M). We define 

74,±(A, M) = Ho(A, M) n g(A, M) x . 

Then 

H(A, M) = Q(A, M) © W ,±(A, M) © H CX c(A, M) . 

We now state the main lemma of this section, which is the key to our 
theorem 

Lemma 4.1. J/-0 e H(A,M) and ij; _L C/(A, M), i/ien /or some ^' e 
H 0t± (A, M) and tp" e H exc (A, M), we have 

P a ^j = */>' + V>" • 
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Moreover, 



(p A iP\p A ij) - y ' \ WW) 

where the ratio is interpreted as zero ifip' = 0. Hence 

7 (A,J,Af)>2J(l-A- 1 )fl- sup M^f 

Proof: First note that P A £(A, Af) = Q(A, Af), soif _L £?(A, Af) then 
P A ^ 1 £(A, Af), which proves that P A ^ G Wo,±(A, Af) © W exc (A, M). 

Now suppose PaV' = ip' + Then Hjjp' = 0. Hence, by Theorem [2 



(p a ip\h~ a p a ip) = (r\H~ A r) > a - a- i )(^'io 



So 



','112 



(P A ^|P A V> " v y V ll^ll 2 , 

By Cauchy-Schwarz and the fact that (ijj'\ip") = 0, we have 

W> = (^'IPaV'} = WW) < WW ■ \\Pa^\ 

I.e. 



< 



ll^ll - ll^ll 

Thus we may define 

*(A,J,Af) = sup 



The positivity of a gap 7(A, J, Af) is then equivalent to 5(A, J, Af) < 1. 
We will numerically estimate S(A, J, Af) in Section 6, directly from its 
definition. For now, we use the lemma to prove the main theorem. 



5. Proof of Theorem 2.2 



We begin this section with the trivial part of the theorem, namely 
the calculation of the spectral gap above the state u] which is the 
translation invariant all up spin state. This is a well-known result, but 
we include it for completeness. To prove that there is a spectral gap, 
we have to prove that there is a number 7j > such that for any local 
observable 

u](5(Xy5(5(X))- l3 5(X*)5(X)) > 0, 



.1 

Ising ' 
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where S(X) = limA/oof-^A, X]. If X e Aa is local, we can take 
8{8{X)) = [#a+[-2 2]> [-^a+[-i i]> ^]] ana - so on - We observe that the 
boundary terms of if J are equal to zero, i.e. 

lim cul(8(Xy(Si L -Sl)8(X)) = 0. 

L— >oo 

So we may rewrite the Hamiltonian in the GNS space of or 
H J = A-'HL + (1 - A- 1 )^ 

oo 

x=— oo 
oo 

-^Ising = (J 2 — Sx^x+l) ■ 

x=—oo 

Clearly H? so > 0, so 7j is bounded below by the spectral gap of H^ sing . It 
is important that u>] is a ground state of both H? so and H^.. It is easy 
to see that the first excitations for Hj si are the one magnon states. 

The one magnon states are those obtained from u] by conjugating with 
observables of the form 

where {c x } is any complex, square-summable sequence. Then one ob- 
serves that u](X) = 0, u](X*X) = 1, and u](X*Hl ing X) = 23. We 
claim that it is easy to see that among all quasilocal perturbations, 
satisfying the first two equalities, these minimimize the Ising energy. 
So the spectral gap of Hj si is 23, and the spectral gap of H J is at least 
(1- A" 1 ) 2 J. It is also very well known that the Heisenberg model acts 
as the discrete Laplacian on the space of one magnon states. We can 
see this since 

rrrj y] — \^ ~ Ca; -l) ^x-1 + ( C x ~ C x+l) S x+1 

L" isotropic ' ^ J / 4 r^p: ^ x i 

x&L 

which implies 

oo 

u>](X* H? BOtIopic X) = J 22 Cx(2c x - c x+ i - Cx-i) . 

x=—oo 

We can choose a sequence of one magnon excitations 
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such that u](XlH^ sotmpic X L ) = 2J/L. Therefore, 

2J(1 - A^ 1 ) < 7j < lu](X* l H j X l ) = 2J(1 - A" 1 ) + 2JA~ 1 L~ 1 , 

for all L, which shows that 7j = 2J(1 — A -1 ). 

The other trivial facts in the theorem are that for the gaps above 
the kink 7(J, M) = 7(J, M + 2J), which follows by translational sym- 
metry, and 7(J,M) = 7(J, — M), which follows by spin-flip reflection 
symmetry. 

We now begin the proof of the nontrivial parts of the theorem. Fix 
J G §N and A -1 in the range < A -1 < 1. 

Definition 5.1. We say that a sequence of triples (L k , M k ,ip k ) satis- 
fies hypothesis (HI) if for all k G N the following holds: L k G N>2, 
M k G [-JL k ,JL k ], i) k in Wo, ± ([l,L fc ] x [l,2J\,M k ), and ||^|| 2 ="l. 
We note that the most important part of (HI) is that ip k G 7io,_i_- In 
particular, this means that P[\,L k \'4 , k 7^ ipk- We say that a sequence of 
triples (L k , M k ,ip k ) satisfies hypothesis (H2) if additionally 

hm ||P[i,L fc ]^fc|| = 1- 

K— >0O 

The main component of our proof is the following 

Proposition 5.2. No sequence satisfies both hypotheses (HI) and (H2). 

We observe of both (HI) and (H2) that if any sequence (L k , M k , ipk) 
satisfies (HI) or (H2), then every subsequence does as well. Using 
this fact and the previous proposition, one can deduce that there is a 
constant 5 > 0, depending on A -1 and J, such that for any sequence 
(L k , M k ,ip k ) satisfying (HI), one has 

limsup ||P[i,L fc ]^fc|| <1-S. 

k— >oo 



Hence by Lemma 4.1 and the discussion following it, 



inf inf 7 ([1, Ll J, M) > 2J(1 - A" 1 )^ . 

LeN> 2 -JL<M<3L 

All the finite volume spectral gaps have a uniform lower bound. Since 
the kink ground states are frustration free, this gives the following 
corollary, which is a reformulation of our main theorem 

Corollary 5.3. For < A -1 < 1 and any J G there is a nonvan- 
ishing spectral gap above all of the infinite volume kink states. 



Proof: (of Corollary |5.3| given Proposition |5.2| ) Let be the 



ground state of the kink. To prove that there is a spectral gap, we have 
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to prove that there is a number 7(J, M) > such that for any local 
observable 

cviU5(xy5(5(X)) - 7 (J, M)5(X*)5(X)) > , 



where S(X) = limA/oof^A; X]. On the other hand, by Lemma |3.4j , for 
any local observable X 

u JM (X) - hm ||^ o(Al)J)M) || 2 > 

where A £ = \—L + 1, L]. If X is a local observable with support in A 
then one can take 6(6(X)) = [H Al , [H Al ,X]], for any A L D A+ [-2,2]. 
Also for M' e [-J|A L |, J|A L |], 

(* (Al, J, M')\(H Al X - XH Al )V (A l , J, M)) = 

because \l/ (A i , J, M') and \l/ (A i , J, M) have the same energy Hence 
5(X)^o(Al, J, M) is orthogonal to all ground states. Thus, we see that 
in this case the infinite volume gap corresponds to the naive guess, i.e. 
the liminf of all finite volume gaps. To be more explicit, let 7(M, J) be a 
positive lower bound on all the finite volume gaps ^([—L + 1, L], J, M), 
which exists by Proposition |5.2j . Define and uj' l to be the states 
on A Al corresponding to \&o(Al, J, M) and 6(X)^ (Al, J, M), respec- 
tively. Then 

^ M (S(X)*S(5(X)) - 7 (J, M)5{X*)5{X)) 

= lim w L (5(X)*5(5(X)) - 7 (J,M)5(X*)5(X)) 

L— >oo 

= lim cu' L (H AL -j(J,M))u; L (5(Xy5(X)) 

L—>oo 

l^oo ||w (A L , J, M)\\ 2 

> 0. 



We will show that (HI) and (H2) are incompatible in a series of 
lemmas. 

Lemma 5.4. If (Lk, Mk,iftk) satisfies (HI) and (H2), then Lf. — > oo. 

Proof: (of Lemma |5.4|) Suppose not. Then there is some sequence 
satisfying (HI) and (H2) and such that Lk = L < oo for all k. But 
Ho([l,L], J) = © A / = _ JL 7^o([l? L], J, M) is a finite dimensional space. 
The finite matrix obtained by restricting and projecting Pun to this 
space is Hermitian, and its largest eigenvalue is 1. Moreover, the 
eigenspace corresponding to 1 is G([l, L], J) = ® A / = _ JL £7([1, £], J, Af). 
All the vectors ipk are orthogonal to £?([!, L], J) by hypothesis (HI). 
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Since finite matrices have discrete spectra, this contradicts hypothesis 
(H2). ■ 

Lemma 5.5. If(Lk,Mk,ipk) satisfies (HI) and (H2) then JLk — \Mk\ —* 
oo. 

Before giving the proof of this easy lemma, we need to define some 
new notation. 

Definition 5.6. For L G N, let Ai(L, J, M) be the set of all vectors 
m = (mi, . . . , m 2 j) G [-L/2, L/2] 2J whose sum is M, i.e. m j — M. 
For each m G M(L, J, M), define 

2.1 

*o([l,^] x [1,2 J\,m) = (g)* ([l,i] x {j}, 1/2, m,). 

Related to this, let Af(J,N) be the set all vectors n = (m, . . . ,ri2j) G 
N 2J satisfying J2j n j = N - Let e = (1, 1, - • - , 1) G N 2J . Recall a 
previous definition 

L], 1/2, n) = "' S r: ■ ■ ■ ■ 

l<x\<X2<-<x n <L 

Define 

2.1 

*' ([1,L] x [l,2J],n) = (g)tf' ([l,L] x {j},l/2, nj ) 

3=1 

Let 7i(N x [1, 2 J], up) be the Guichardet Hilbert space 

(g) (C 2 ^, 1+1/2)^)). 

(x,j)6N[l,2J| 

Define V([1,L] x [1,2J],N) to be the projection on H(N x [l,2J],up) 
which projects onto vectors ip such that 

(x,j)e[l,L]x[l,2J\ 

and Sf xJ) ij = for (x,j) G (N\ [1,L]) x [1,2J\. Finally, let 

%([1,L] x [1,2 J], N) 

be the unique normalized ground state of -f/n l]x[i,2J] ira range of 
V([1,L] x [1,2J|,JV). Specif toe phase to that %([1,L] x [1,2J\,N) 
has real coefficients in the Ising basis. This ground state exists and 
is unique (and has real coefficients) since Hun x uzl acting on the 
range £>([!, L] x [1,2J\,N) is unitarily equivalent to the finite matrix 
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l]x[i,2J] acting on 7Y([1, L] x [1, 2 J], M). (We apologize for the abuse 
of notation : the same notation is used for the operator acting on two 
different, but isomorphic, Hilbert spaces.) 

The following are easy and useful observations. The set 

{*o([l, L] x [1, 2 J], m):meM(L, J, M)} 

is an orthogonal basis for H, ([l,L] x [1,2J],M). The set of indices 
J\f(3, n) is finite, while if one defined the analogue of A4(L,3,M) re- 
placing L by N it would not be finite. There is a simple translation 
between the two vectors defined above: 

%([l,L] x [l,2J],n) j^zMMj x [MJjjLe-n) 
\\%([1,L] x [l,2J],n)|| ~ ||* (M] x [l,23],\Le-n)\\ 
By Lemma |3.4| , the following strong limit exists 

*' (N x [1, 2J], n) = lim *' ([1, L] x [1, 23], n) . 

L— >oo 

One has the following simple formula for the action of P([1,L] x 
[1,23], N) on H ([l,L f ] x [1,2J],M): 

V([l, L] x [1, 2J], iV)*' ([l, L'] x [1, 2J], n) = %([1, L' A L] x [1, 2J], n) . 
By Lemma |3.4| , again, the following limit holds for all n £ A/"(J, iV): 



V 7 £-ooL'-»oc ||*' ([1,L'] x [l,2J],n)|| 

By Lemma |3.4|, again, the following strong limit exists 



(5.5) %(N x [l,2J],iV) = lim %([l,L] x [1,23], N). 
We note for the reader that for large enough L 

%([1, L] x [1, 23], N) ex *o([l, ^] x [1, 2J], n) , 

neJ\T{3,N) 

and that 

$i(Nx [l,2J],iV)cx ^ f (Nx [l,2J],n). 

neJ\T{3,N) 

The proportionality constants are necessary because both \l/g([l,L] x 
[1, 23], N) and ^(N x [1, 23], N) are chosen to be normalized. 
Proof: (Lemma |5.5| ) We will verify that — Mk —> oo; the fact 
that 3Lk + Mk — > oo then follows by spin-flip/reflection symmetry. If 
it fails for some sequence, then some subsequence, also satisfying (HI) 
and (H2), has the property that 3Li k — M\ k = M is a finite constant. 
We assume that this subsequence was taken at the beginning to avoid 
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double subscripts. Since each ipk G 7~£o,±([l> Lk] x [1,2J],M), we can 
write 

V- t \ *o([l,^fc] x [l,2J],|L fc e-n) 



||*o([l,^fc] x [1,2J], |L fc e-n) 



So 

*{,([!, x [l,2J],n) 



ntMVM.m J*i([l,iJx[l,2J|,„)ir 

We know L k — > oo by the last lemma. So ^([l,!^] x [l,2J],n) — > 
^(N x [l,2J],n) as A; — > oo for each n. Furthermore, {cfc(n) : n G 
M(J,N)} is a unit vector in the finite-dimensional space CUWOI, for 
each fc. So there is a subsequence with Q fc (n) — > c(n) for each n. 
Again, we assume this subsequence was chosen at the beginning to 
avoid double subscripts. Thus, ip' k converges to the vector 

tf,= c(n)*' (Nx[l,2J],n). 

For L < oo and k large enough, P[i,l] > ]■ Hence, by (H2) 
(ip\P[i,L]ip) = lim {il)' k \P {lt L}ij)' k ) > lim (^l^i,^]^) = 1 

for all L. Also, £T[i,L fc ]x[i,2J] > #[i,£]x[i,2j], so #[i,l]x[i,2J]?/' = 0. Recall 
~ i 

P[l,L]H[i t L}x[l,2J}P[l,L] = -^jP[l,L]H[l,L]x[l,2J]P[l,L] ■ 

So 

(5.6) ■H"[i,i,]x[i ) 2J]V' = 

for every finite L. 
Notice, 

(5.7) 

lim ||X>([1,L] x [l,2J],iV)V>|| = lim lim ||X>([1,L] x [1, 2J], N)ip k \\ = 1 

L— >oo L— >oo fc— >oo 



by (Ml ). Together, (gj) and (p?j\) imply that 

lim |(^|^([l,L]x[l,2J],iV))| = l 

i.e. that ^ = e^^N x [1, 2J], N). Since ^ -> 
(5.8) h lB*W = , 
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On the other hand, we know (ip' k \^' ([l, L k ] x [l,2J],iV)) = 0, because 
tp k JL G([l,L k ] x [1,2J],JV). So, by this implies 

(5 . 9) Bm IWW(Nx[l,2Jl,JV))| =0- 



Clearly, ( |5.8|) and (5.9) are incompatible and we have a contradiction. 



Proof: (Proposition |5.2|) 



Let .M (J, M) be the set of all (m b . . . , m 2 j) with . rrij = M. De- 
fine for any m G Z, 

^([-a,6],l/2,m) = ^ £ 

n — o<ai<— <a n <0</3i<— <;9 m+ „<6 
/3 1 + .../3 m+n -a 1 -...-a n>s + . . . g+ g ~ . . S ~ Q z . 

Let 

2J 

*(,([-a,6] x [l,2J],m) = (gK([-a,6] x {j}, 1/2, m,) . 

i=i 

Suppose that (L fc , M k ,ip k ) is a sequence satisfying (HI) and (H2). By 
Lemmas [5.4| and |575| we may assume that L k — > oo and JL k — \M k \ — ► oo. 
We may choose a subsequence so that JL^ — = 2Jr^ — A 7 for all k, 
where r k G N>o is arbitrary and N G [0, 2 J — 1] is fixed. Let 

let Afe = [1 — r k , L k — r k ]. Then we can write 

%(A k x [l,2J],m) 



\%(A k x fl,2Jl,m)|| ' 

Now, for each fc, {Cfc(m) : m 6 M(J, N)} is a normalized Z 2 sequence, 
instead of a finite dimensional vector. So we only know that a weakly 
convergent subsequence exists, not a strongly convergent one. We need 
some kind of tightness result, which is given in the following lemma. 

Lemma 5.7. Let Mr(J, N) be those m G M(J, N) such that all com- 
ponents lie in the range [-R + 1, R]. Then we have 

lim liminf > |Cfc(m)| 2 = 1. 

R— >oo fc^oo ' J 

meM R (J,N) 

We give the proof of this technical but important lemma at the end 



of the section. First, we see how Proposition E72 follows. 
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Lemma 5.8. Define Tr to be the "projection onto those vectors with 
all down spins at sites (a,j) when a < —R, all up spins at sites (a,j) 
with a > R + 1, and Sf_ R+1 fl i x r 1 2 j\ equal to N. We have 

lim liminf ||.7 7 rV4|| = 1 . 

R—*oo fc— »oo 



Proof: By lemma for any e > 0, we can choose R large enough 
that 

liminf |C fc (m)| 2 > 1 - e . 

k. — >r>r> * * 



k-^oo 

m£Mn(J,M) 



By Lemma |3.4|, the following strong limit exists 



%(Z x [l,2J],ra) = lim %(A k x [l,2J],m) 

fc^oo 



Furthermore 



||.T fll tt' (Zx[l,2J],m)|| = 
Ri-oo ||$' (Z x [l,2J],m)|| 

Thus 

liminf H^Vfell 2 > 

-Ki— >oo 

Since e was arbitrary, we are done. I 

We know that H-^ijj' k = and lim^ H-Pa^^H = 1. As before, these 
two facts are sufficient to guarantee that for any finite A C Z, 

hm%i = 0. 

fc-+oo 

In other words, letting G\ be the projection onto the ground state 
space of Hj^, that 

(5.10) lim ||G A ^|| = 1. 

fc— >oo 

We note that for A = [-R+1, R], the projections Tr and Ga commute, 
and in fact their product is the projection onto the normalized ground 
state vector 

%([-R+l,R] x [l,2J],iV). 
By Lemma |5.8| and equation ( |5.10| ) we see that 

lim liminf \(%([-R+l,R] x [1, 2 J], N)\i/;' k )\ 2 = 1 . 

K—>oo fc— >oo 

Since V' ([-R+l, R]x[l, 2 J] , N) converges in norm to %(Zx[l, 2 J] , N) 
as R — * oo, we have 

(5.11) liminf | (^(Zx [l,2J],iV)|^)| 2 = 1 . 

fc^oo 

Now comes the contradiction. We know, by virtue of the fact that 
i// k e H 0i± ([l,L k ] x [l,2J],M fc ), that f k _L %([l,L k ] x [l,2J],iV). 
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But on the other hand, we know that ^([1, Lk] x [1, 2J], N) converges 
strongly to %(Z x [1, 2J], N), so 

limsup|(^(Zx [1,2J],JVM>| 2 = 0, 

fc— >oo 

clearly contradicting ( [5.111 ). | 



Proof: (of Lemma |5.7|) We define -Dr to be the projection onto all 
those vectors in 7i(Z x [1,2 J], kink) with at most 2J — 1 down spins 
shared between the sites (-R+ 1,1),..., (-R+ 1, 2J), and at least one 
down spin shared between the sites (R, 1), . . . , (R, 2J). Let 

*o(A,m) " 



ll«o(A,m)|| 
Let A R = [-R + 1,R\. It is clear that if 

i)= C(m)y (A R ,m) 

m£M(J,M) 

then 

(^| J D^) = ^|C(m)| 2 || J D^ (AR,m)|| 2 

m 

> Yl \C(m)\ 2 \\D R MA R ,m)\\ 2 

meM(],M)\M R (3,M) 

> \ E i^)i 2 

meM(3,M)\M R {J,M) 

the second line owing to the fact that D R commutes with S?_ L+1 n x j,} 
for each leg j and all L G N. The last inequality is the result from 
to the following consideration: if one of the rrij is greater than R or 
less than — R + 1, then with probability at least 1/2 one finds the state 
with at least one down spin at R, or at least one up spin at —R + 1, 
respectively. Thus 

£ |C(m)| 2 > 1-2^\D R ^). 

meM R (3,M) 

Now suppose that ip is a ground state vector of ifj . Then, because 
D R commutes with S~ we have 

Oj\D r ^} < max {V (A R ,M)\D R y (A R ,M)). 

M£[-2JR,2JR] 
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We can make the following crude but simple estimate 

~ 4 T 2 ~Ra 2R 

(5.12) max ||Z>A(A fl , M)|| 2 < 



Me[-2JR,2JR] ' 1 - AJ 2 Rq 2R ' 

We prove this estimate for M < 0, and the M > follows by spin 
flip/reflection symmetry. Of course if M = —2JR, then the estimate 
is trivial because \l/o has all down spins, so the expectation with Dr is 
zero. Suppose — 2JR < M < 0. Then there is at least one down spin 
and at most 2JR down spins. We write M(A#, n) for all the classical 
Ising configurations on the lattice Ar with exactly n down spins. We 
write M(A^,n, j) for all the configurations with the extra constraint 
that the number of down spins on the sites (R, 1), . . . , (R, 2 J) is j. The 
normalized ground state \I/o(A_r, 2JR — n) is 

*o(A B ,2Ji? -n) = Z- 1 Yl W({m(a,j)})\{m(a, j)}) , 

{ro(a,j)}eM(A Bl n) 



where 



and 



W({m(a,j)}) = ] I q~ am ^ 

(«>i) 



Z 2 = E W({m(a,j)}y 

{m(a,j)}eM(A R ,n) 

On the other hand, defining Erj to be the projection onto those states 
with exactly j down spins on the sites (R, 1), . . . , (R, 2J), we have 

Z 2 

\\ErjV ([-R + 1,R] x [l,2J],2Ji2-n)|| 2 = -^ 

where Zj is the same as Z but with the sum over classical configurations 
restricted to ~M(A R ,n, j). 

The inequality comes from recognizing that Z 2 +1 < 4J 2 Rq 2R Z 2 . This 

is a straightforward estimate. Define a lexicographic order on Ar by 
(oil, ji) < ( a 2,j2) if either ot\ < a 2 or a.\ = a 2 and j\ < j 2 . Define 
a map / : M(A_r, n, j + 1) — > M.(Ar, n, j) where the down spin at 
the greatest site (a, j) with a down spin, is exchanged for the up spin 
at the least site (/?, k) with an up spin. Note that since there are at 
most 2JR down spins, the point (j3, k) must lie in the subset [-R + 
1,0] x [1,2J], so there are at most 2JR choices of (/3,k). Similarly, 
since j + 1 > 1, we know there is at least one down spin in the sites 
(R, 1), . . . , (-R, 2J). So there are at most 2J choices for (a, j) = (R,j). 
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Thus, #/~ 1 ({m(a, j)}) < 4J 2 i? for any configuration {m(a,j)}. We 
see that 

W(f({m(x,j)})) > q- R W({m(x,j)}) 

because the down spin at site (a, j) = (R, j) has moved at least R units 
to the left to (f3,k), /3 < 0. So 

Z 2 1+1 = E W({m(a,j)}Y 

{m(a,j)}£M(A a ,n,j+l) 

<q 2R E ^(/({m(a,j)})) 2 

{m(«J)}eM(A fl ,nj+i) 

<4J 2 i?g 2i? E ^({m(a,j)}) 2 

{m(a,j)}6M(A fl ,nj) 

= 4J 2 Rq 2R Z 2 . 

The crude estimate is proved. 

Now we know that as k — > oo, the vectors come closer and closer 
to the ground state space of . So our estimate implies 



liminf V \C k {m)\ z > 1 

mG.M fl (J,M) 

Hence 



8J 2 i?g 2ii 

fcimir i ^ - 

1 - AJ 2 Rq 2R 



2 D„2ii 



hm hmmf \Ch(m)\ > 1 — hm ^ 

R^oo k^oo ^ ' V n ~ R^oo 1 - A3 Ra 2R 
m&M R {J,M) 

and this concludes the proof. 



6. Numerical Approximation 

We now find an explicit representation of 

Proj(ft (A,M))P A Proj(ft (A,M)). 

From this we numerically calculate 1 — S(L, J, M). We begin with some 
definitions. First of all, we will always have A = [1,L] in this section, 
and hence A = [1, L] x [1, 2J]. For N G [0, 2JL], define 

P(K, N) = Proj(H (A, JL - iV))P A Proj(^ (A, JL - N)) . 
Also, define the "classical Ising configurations" to be 

M(L,2J,A0 = {AG[0,1] X : ^ A(x,j) = N}. 

(x,j)eA 
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These are {0, l}-matrices with 2 J rows, L columns, and N ones. For 
any A G M(L, 2 J, N) we define 

(x,j)eA 

where Q = <S)( X j)eA I V2)o,j)- Then it is clear that TC(A, JL — N) has 
an orthonormal basis {<pA '■ A G M(L, 2J, N)}. For any matrix A, we 
define two vectors r a E C 2J and c A E C L by 

L 2J 

r A (j) =^2A(xJ) , c A (x) =^2A(xJ) . 

x=l j=l 

Finally we define 

(6.13) M r , c = #{A G M(L, 2J, iV) : r A = r , c A = c} . 

Note that by its definition M rjC is unchanged if one permutes the com- 
ponents of r or c. We mention that there is no known formula for 
M rjC although it has useful characterizations in terms of generating 
functions. (C.f. j27| §7.4 for more details.) 

The definitions immediately lead to the following result. 

Lemma 6.1. The following are true identities: 

(6.14) vP (A, J, \Le - n) = ]T <j> A <r aA 

AeU(L,2J,N) 
r A =n 

(6.15) ||* (A, J^Le-n)\\ 2 = £ q 2xc M n , c 

ce[o,2J] L 



\ \ca(x) 



x=l x " v " BeU(L,2J,N) 

C-B=C-A 



(6.17) 



(^ (A, J, ^Le - m)|P A ^ (A, J, ^Le - n)) 



' 2J 



£ M m , c M n , c q 2xc f[f 

lt).2.J\ L x=l ^ 



ce[o,2j] 



c(x) 
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We can define an action of ©2j on H,q(A, JL — N) by U {tt)-4>a = 4>kA, 
where 7rA(x,j) = A(x, n~ l (j)). By fl6.17| ), we know that the range of 
P(A,N) is a trivial representation of ©2j- We define F (L,2J, N) to 
be the set of all sequences \i = (/ii, . . . , /i 2 j) such that 

L > /ii > Hi > ■ ■ ■ > fi2j > 

and fii + ■ ■ ■ + H2j = N. These are restricted partitions, but allowing 
parts equal to zero. (We mention this fact for consistency. The inter- 
ested reader can consult 0, fTTj| , or |^7| for more information about 
partitions.) Then the range of P(A, N) is spanned by the linearly in- 
dependent vectors 

^ = T^j E U (*)*(A, XL - //) , fJ, G P (L, 2J, N) . 

^ 7TS62J 

We note that we can define an action of ©2j on Z 2J in the obvious way, 
with the outcome that 

\P(A, JL - vrn) = C/(tt)^(A, JL-n). 

The orthogonal basis {ip^ : fi G Wq(L, 2 J, A)} is not yet orthonormal. 
We observe that 

||^|| 2 = ((2J)!)- 2 ||v[/(A,JL- / i)|| 2 

x # (orbit of jj) x # (stabilizer of /i) 2 . 

Since # (stabilizer) x # (orbit) = #©2j = (2J)!, we have 
(6.18) 

||^|| 2 = #(orbit of ii)-i\\${k,ZL-ii)f 

f 2J 
\n (ii), ni(/i), n 2 (/i), ...,n L (ii) 

E«c(*)=Jv 

where rik(fi) is the number of parts of /i equal to fc. 

By ( |6.17|) , and the fact that the inner product is invariant under 
permutations of m and n, we obtain 

L 

{%\PM= M m , c M niC q 2x - c l[ 

ce[0,2J] L x=i 
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for all //, v G Po(£, 2J, A r ). Therefore P(A,N) is represented, on its 
range, by a matrix = (V(fi, v) : fj,, v G Po(£, 2J, A r )), where 

2 J \ 1/2 / 2 J \ 1/2 

,n (/f), • • • , n L (/i)y \ n o(^), ■ ■ ■ , ^l(^)/ 
/ \ 



x 



(6.19) 



:ffl 9..Tl i 1=1 ^ ' 



cG[0,2J] 



/ 



E 



cG[0,2J] L 
\£ x c(x)=N 



J 



This is quite a complicated looking formula. There are two nice fea- 
tures about it. First, the powers of q grow quite rapidly. Second, the 
matrix [M^y : /i, v G Po(£, 2 J, JV)) is upper triangular with respect to 
dominance order, where v' is the transpose of the partition v. (Domi- 
nance order is the natural partial order on partitions.) These two facts 
insure that the matrix components of V decay exponentially with the 
distance from the diagonal. (More details and an equivalent expression 
are found in p8| .) 

We have not attempted a rigorous analysis of V, but we have ob- 
tained very convincing numerical data, by simply numerically diagonal- 
izing the matrix. The main qualitative feature of the lower bound for 
the spectral gap is that it is not always maximized at the Ising limit. 
In particular, if J > 3/2 and the number of down spins N satisfies 
N = 2J [|A|/2j, then the local maximum for the lower bound of the 
spectral gap occurs somewhere other than the Ising limit. The Ising 
limit is a classical model, whose energy levels can be calculated explic- 
itly, and doing so it is easy to see that the lower bound for the spectral 
gap is actually equal to the true spectral gap at the Ising limit. So the 
true spectral gap of the XXZ spin chain with J > 3/2 has a local maxi- 
mum somewhere other than the Ising limit in finite volumes as long as 
the number of down spins satisfies iV w J|A| and N = 0(mod 2J). One 
can take N ~ J|A| and N = 0(mod 2 J) because of the approximate 
periodicity of the spectral gap in N. This same result is also obtained 
by Ising perturbation series, where we show that the curve for the spec- 
tral gap is concave up at the Ising limit. Also, the asymptotic analysis 
of Section 7 verifies the qualitative picture of the spectral gap when 
J > 1. 
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We now give an alternative description of V in terms of representa- 
tions. Recall that the ground state space of the spin 1/2 XXZ model is 
the highest dimensional irreducible representation of SU 9 (2). On the 
other hand, the symmetric tensors form the highest dimensional irre- 
ducible representation of SU(2). Consider the array A = [1, L\ x [1, 2J]. 
At each site put a two dimensional representation of SU(2) and SU g (2). 
Note that this is possible because the two dimensional representations 
of SU(2) and SU g (2) coincide. Now tensor all the representations in 
a single row, considering them as representations of SU g (2). For each 
row, define an operator Rj(q) which projects onto the highest dimen- 
sional irreducible representation of SU 9 (2). Next tensor all the rep- 
resentations in a single column, considering them as representations 
of SU(2). Define C x to project onto the highest dimensional irre- 
ducible representation. Then the operator ^PfA, N) is identical 
t° rL=i C x llj=i ^7(9); modulo null spaces. If one turns the procedure 
around, first projecting on columns then on rows, one almost (but not 
quite) recovers the original problem of the spin J XXZ chain. 

6.1. Perturbation Series about Ising Limit. We now perform a 
perturbation analysis for 7QI, L], J, M, A" 1 ) about the point A^ 1 = 0, 
i.e. the Ising limit. We write 

H(A~ l ) = H<® + A- 1 ^ 
tf (0) = £(J + ^)(J-^ +1 ) 

x=l 

L ^ ( 1 _ 1 

x=l ^ 

We have left out of the first order corrections to the boundary 
terms. However since all our vectors are local perturbations of an 
Ising kink, the first order corrections to the boundary terms will act 
as a multiple of the identity. We will include these trivial corrections 
after we perform the perturbation theory with as above. We note 
that ^(A" 1 ) is unitarily equivalent to H{— A -1 ), where the unitary 
transformation is 

/ \L/2\ 

U = exp 2m 5 2 3 

V j=i 

This proves that the point A^ 1 = is always either a local maximum 
or a local minimum of 7QI, L],J, M, A -1 ). If 7QI, L], J, M, A" 1 ) is 
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differentiable near A -1 = 0, the first derivative is zero, and we proceed 
to second order perturbation theory. The reason 7QI, L], J, M, A -1 ) 
may not be differentiable near A -1 = is that the first excited state 
may be infinitely degenerate in the Ising limit. This is the case for spin 
1/2 and for J = 1 when M is odd. For J = 1/2 and any igZ there is 
an Ising ground state 



*o(a0 



It is easy to see that there are infinite families of first excitations, for 
example Ylf =1 S^ +1 _jS~ + j_{ i So(x) for any L > 1 and y > x + 2 - L. 
For J = 1 and M odd, the ground state is 

^y<x / \y>x 

and there are two classes of excitations, each infinitely degenerate: 
SyS~^fo(x) for any y < x; and S^S~^o(x) for any y > x. These 
are the only cases where the first excitations are infinitely degenerate. 
The only cases where there is a finite degeneracy for the first excited 
state are J = 2,3,4, .. . and M congruent to J modulo 2J. Then the 
ground state is 



*o(a0 



and the two first excitations are S^_ 1 S~^/ (x) and S£ S~ +1 ^J (x). The 
other most interesting case, which has unique first excitations are J > 1 
and M divisible by 2 J. Then 

*o(x) = 

and the unique first excitation is S£ S~ +1 ^o(%)- There are ground states 
with infinitely degenerate second excitations, and so on, but this does 
not interest us. 

We now consider the results of second order perturbation theory, 
assuming the first order excitation is non-degenerate. The kink ground 
states of the Ising model are all of the form 



Mx,n) = I Q9|-J)„ ) ® 

,y<x 
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where one can assume that < n < |_JJ • The first excited state is then 




Vi(x,n) = ( 05 | -J) j, I ® |-J + n + l) x ® |J-l)x+ 
,y<2 / 

which has energy = n + 1. We now expand to determine the 
corrections for small but nonzero A" 1 . In particular we write E(A~ 1 ) = 
E(°) + A^E^ + .... The perturbation series is standard, so we omit 
details. The results are that E^ = and 

This is not an accurate description of the kink Hamiltonian because 
we have not included the correct boundary fields. To fix this situation 
we must add J(\/I — A~ 2 — l)(Sf — S^). It is obvious that for a long 
enough spin chain, and excitations which are localized at the interface, 
the extra boundary fields act just as — 2J 2 (Vl — A~ 2 — 1) times the 
identity. Note that this is A" 2 J 2 + o(A~ 2 ). So for < M < J, 

7 ([1,L],J,M,A- 1 ) 

A-1=0 



(6.20) 



M t (J+1)(2J-1) 2J 2 



2 M + 3 2J - M - 1 

The finitely degenerate case is J G Z and M = J mod (2J), as men- 
tioned before. Then the first excitation of the Ising ground state is 
doubly degenerate, and we perform degenerate perturbation theory. 
As soon as A~ x > 0, the degeneracy lifts and there are two branches. 
It is easily verified that the curvature of both branches is negative, but 
we are only concerned with the lowest branch which gives 



d 2 

(6.21; 



7 ([l,L],J,J,A- 1 ) = -8--^- 



J-l 2 J+3 



We list some values for the curvature of 7 in Table |6.1| . In the Ising 
limit, the minimum gap occurs for n — 0. One can see from this table 
that for M = and J > 1, the gap is concave up at A^ 1 = 0. This is 
the basis for Conjecture 



7. Boson Model 

We now give a heuristic derivation of the free Bose gas model for the 
XXZ spin system in the limit J — > 00. This is an approximation to the 
full XXZ Hamiltonian on a finite chain A = [1, L]. Our approach 



is similar to that of [13], although we would suggest to the reader to 



look at H 0] and [|| instead. We are interested in the classical limit, 
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n = 


1 


2 


3 


J= 1/2 


— oo 








1 


-1/3 


— oo 






3/2 


11/12 


-9/4 






2 


7/3 


-1/4 


-46/5 




5/2 


97/24 


4/3 


-39/20 




3 


91/15 


3 





-26/3 



Table 1 . Some values of the curvature of the gap in the 
Ising limit. For J = 1/2, n = 0,and J = 1, n = 1, the 
excited state is infinitely degenerate, and the curvature 
is infinite as well. 



J — > oo. Lieb, ]nj, proved that for the Heisenberg model one obtains 
the classical partition function as a scaled limit of quantum partition 
functions. Lieb's method used coherent states to obtain rigorous upper 
and lower bounds on the partition function. In || , the same result was 
derived without coherent states, and then it was shown that with a 
sufficiently large external magnetic field the large J limit of the XXX 
model can be viewed as a free Bose gas, which verified predictions of 
Dyson in |], |10[. A more recent proof of the Bose gas limit has been 
obtained by Michoel and Verbeure [E3], where the spin J operator is 
viewed as a sum of 2 J spin 1/2 operators, (using a spin ladder), and 
then a noncommutative central limit theorem is applied. 

The physical requirement of a large external field to obtain the Bose 
gas limit in the isotropic case is easy to understand. All the approxima- 
tions (including our own) rely on a spin wave description of the elemen- 
tary excitations. In order for this to be valid, the ground state must be 
very nearly saturated, i.e. the ground state should satisfy | (S^) | — J <C J- 
To accomplish this for the isotropic Heisenberg model, one must place 
a rather large external magnetic field. This is an important difference 
between the isotropic and anisotropic ferromagnets. The XXZ model 
with A > 1 possesses kink ground states, which the isotropic model 
does not. For the kink ground states, there is a quantum interface sepa- 
rating two regions, which we can assume to be located at a = 1/2. For 
sites a < one has (S^) < -J + Cq~ a and for a > 1, (S*) > +3-Cq a . 
Thus the spin is saturated well away from the interface, with exponen- 
tially small corrections. The XXZ model exhibits saturation with just 
a boundary field, and the boundary field vanishes in the thermody- 
namic limit. Moreover the boundary field is known to give the correct 
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ground states (cf |16|1 ). So the boson picture is quite natural for the 
XXZ model. 

The energy-momentum dispersion relation is different for the XXZ 
model than for the isotropic model, as one would expect. The most 
important difference is that the lowest energy spin wave is not actually 
localized in momentum space, but in position space. It is localized 
at the interface, instead of being spread out uniformly over a large 
region. (There is one other spin wave with lower energy, in fact zero 
energy. But this is the spin wave which simply moves one ground state 
to the other, owing to the fact that all ground states in all sectors of 
total S 3 have equal energy. We remove this boson by restricting to 
a single sector.) Moreover, there is a spectral gap between the lowest 
spin wave, and the others. The next independent spin wave boson does 
have a well-defined momentum, and from there on the usual picture 
of spin waves prevails. These are the results for one dimension, but 
the Bose gas model also holds for excitations of the (1, . . . , 1) interface 
ground states in dimensions d > 2. In dimensions higher than one, 
the low lying spectrum is more complex, having a continuous band of 
interface excitations at the bottom, as proved in [^], |J. An interesting 
recent result by Caputo and Martinelli gives a rigorous lower bound 
for the spectral gap in a large but finite system A, whose power law is 
| y\ I — j n a g reem ent with ||. This gives strong evidence that the only 
excitations beneath a certain energy are interface excitations. The 
Bose gas approximation implies more, that for large J the interface 
excitations are separated from all other spin wave excitations by a 
spectral gap of order J. The d > 2 results will be elucidated in a 
forthcoming paper [^], as will be the rigorous proof of the spin wave 
Boson model for the XXZ model. For now we provide a heuristic 
argument. 

We begin by considering the simplest case, namely A = [1,2]. It 
is convenient to work in the dual space to 7i(A, J). Namely, for ip : 



[-J, J] x [- J, J] C, define 





mi,rri2G[— J, J] 
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Then Hj^ip) = \Gjip) where Gj is an operator on C' _J,J ] 2 
GjV>(mi,m 2 ) = (J 2 - mim 2 + A(A)J(mi - m 2 ))ifj(mi, m 2 ) 
-±[J(J + l)- mi (m 1 + l)]^ 

[J(J + 1) - m 2 (m 2 - l)] 1/2 ^(mi + 1, m 2 - 1) 

-^[J(J +l)-mi(mi-l)] 1/2 

[J( J + 1) - m 2 (m 2 + l)] 1/2 ^(mi - 1, m 2 + 1) 

where A(A) = s/l — A~ 2 , and we define ip(mi,m 2 ) = for any 
(mi,m 2 ) not in [— J, J] x [— J, J]. All we have done is to explicitly 
write down the action of the spin matrices. Next, for any real num- 
bers — 1 < fii < 1, i = 1,2, let us define a linear operator T W)M2! j : 
C°°(R 2 ) -> Cl- J ' J l 2 , where T^,.,* = 

ijj(m 1 ,m 2 ) = *(J~ 1/2 (m! - mJ), J~ 1/2 {m 2 - fi 2 J)) . 

This operator has a very large null space. But if for fixed 1^1,^2 one 
knows that T MliM2i j\l/ = for every J, then \1/ must obviously also be 
zero. 

There are many choices of operators CKj on C°°(IR 2 ) which satisfy 
GjT^ ltfJi2t j = Tfj, ufX3t jJ{,j. One particularly good choice is the following 

(7.22) 

X 2 ) = J 2 ^1 - (fa + Xi J~ 1/2 )(^2 + ^2-T 1/2 ) 

+ A(A) ({m - (i 2 ) + (xi - x 2 )J- x / 2 ) J*(sci,ac 2 ) 
~2A" S ( II [j(J + l)-(M + ^J 1/2 )(M + ^J 1/2 + (-l) m e)] 1/2 

£=±1 \ 1=1,2 

* (xi +eJ" 1 / 2 ,x 2 -eJ" 1 / 2 ) ^ , 

which is the same as the definition of Gj, but now allowing the operator 
to act on smooth functions instead of discrete functions. One can 
formally expand 

ttj = J 2 X^ + J 3 / 2 <K^ + JftW + . . . 

considering the shift by iJ" 1 ^ 2 as e ±J 1/2,9 , and expanding in the small 
parameter J" 1 / 2 . The resulting expressions for "K^ 2 \ "K^ 3 ^ 2 ^ and "K^ 
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are as follows 
(7.23) 

^ = 1 - Ml /x 2 + A(A)( Ml - /la) - A- X [l - m?] 1/2 [1 - ^] 1/2 
(7.24) 

5{( 3 / 2 ) = A(A)(x a - x 2 ) - Hix 2 - /i 2 Xi 



+ A- 1 



/if 

01^1 H , /^2^2 



V 1 -01 0^02 



(7.25) 



X 2 X 2 . 2/ii/i 2 XiX 2 



+ 



(i-/i?) 2 (i-m 2 )(i-mI) 

1 1 



Let us now use the parameter 77 = log(l/g), which is related to A by 

A -1 = sech(?7) , A(A) = tanh(?7) . 

Then ([7T23D shows that "K^ is a multiplication operator, multiplying 
by the non-negative constant 

sechfo) [fl 2 y/(l + /iO (1 - /i 2 ) - e""/ 2 ^(1 + /x 2 ) (1 - //i)) ' ■ 
Therefore, J{( 2 ^ = iff 

(7.26) 3r e R s.t. Va G A, /t Q = tanh(r/(a - r)) . 

The number r is determined by fii + /t 2) implicitly. Note that IK 1 - 3 / 2 ) 
is also a multiplication operator, but given ( [7.26Q we know that it 
vanishes identically, as well. So the first non- vanishing term is Ji^. 
We use (3) to rewrite (|7T25| ) 

IK (1) = — sech(r/(l - r)) sech(r/(2 - r)) 

x ( - <9 2 + cosh 4 (r/(l - r))x? - <9 2 + cosh 4 (r/(2 - r))xl 

(7.27) V 

+ 2d xl d X2 — 2 cosh 2 (?7(l — r)) cosh 2 (?7(2 — r))xix 2 
— [cosh 2 (?7(l — r)) + cosh 2 (?7(2 — r))l 



Now we notice the following: CK*- 1 ) is a second order differential oper- 
ator which is homogeneous in d Xl , d X2 , X\ and x 2 except for a constant 
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(the zero point energy). Therefore, "K^ can be regarded as the Hamil- 
tonian for a two-mode Boson system with quadratic interaction. Thus, 
for a = 1, 2, we define 

a a = -i= (cosh(r/(a - r))x a + sech(?7(a - r))d Xa ) , 

(7.28) Y 

al = -j= (cosh(?7(a - r))x a - sech(r/(a - r))d Xa ) , 

which satisfy the Canonical Commutation Relations 

(7.29) [a a , a p ] = [a* , ap] = , [a a , a£] = <J ai/J , 
and also 

(7.30) a^a = 2 ( cos h 2 (^(a — — sech 2 (?7(a — r))d Xa — l) . 
One sees that 

a n + al = v / 2cosh(?7(n — r))x n 

and 

a n - al = v / 2sech(?7(n - r))d Xn . 
These relations imply 

2 cosh(?7(l — r)) cosh(7?(2 — r))x\X2 = (Si + a\){a,2 + a^) 

and 

2sech(?7(l -r))sech(?7(2 - r))d Xl d X2 = (Si - S{)(S 2 - aj) > 
respectively. Hence 

29,., <9 X9 — 2 cosh 2 (?7(l — r)) cosh 2 (?7(2 — r))xix 2 

(7.31) 

= —2 cosh(?7(l — r)) cosh(r/(2 — r))(S|S 2 + S3>Si) . 
All of these algebraic manipulations allow us to rewrite ( [7.271 ) as 
(7.32) 



A 



cosh(?7(l — r)) A + A cosh (r/(2 — r)) * ^ ± 
ajai H — a 2 a 2 — a|a2 — a 2 ai 



cosh(?7(2 — r)) cosh(?7(l — r)) 

One can verify that one of the two eigenmodes of this system has zero 
energy, while the other is positive. The zero mode is a direct conse- 
quence of the infinitely many ground states, corresponding to different 
values of M, the third component of the total spin. If one considers the 
canonical picture, restricting the total magnetization to have a fixed 
quantity, then that boson disappears. 

The above results extend directly to a chain of arbitrary length. One 
has jLti, . . . , fi£ satisfying ( |7.26| ). This is required for J-C^ to vanish, and 
sufficient for IH^ 3 / 2 ) to vanish. The definition of the single site Bosons 
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is just as in (|7.28|) , for each a G A = [1, L], and the definition of 3fW 
becomes 

(7.33) DC^ = Jvjalflp, 

where J a ^ = if \a — (3\ > 1, jTa^ = — A -1 if |a — (3\ = 1, and 
_i cosh(?7(l — r)) 



cosh(?7(2 — r)) 

Jl,l = A " 



x cosh(?7(L — r)) 



\Ja,a A ^ 



cosh(?7(L — 1 — r)) 
cosh(?7(a — r)) 



cosh(?7(a + 1 — r)) 



cosh(r?(a; — r)) 
H w / ^tt , for 1 < a < L . 



cosh(?7(a — 1 — r)). 

The single site Bosons are coupled, but the coupling matrix can be 
diagonalized. We obtain new Bose operators b n , in terms of which the 
Hamiltonian is diagonal quadratic: 

L-l 
n=0 

where b n = ^2 a Va ? a a with the eigenvector of J corresponding to 
the eigenvalue X n . Again there is one zero-mode, A = 0, and Aj > 0, 
for 1 < % > L — 1. Ai remains isolated in the limit L — >• oo, for 
< q < 1. 

By considering only the leading order terms for the bottom of the 
spectrum, we have exchanged a quantum many body Hamiltonian Hj^ 
to an L body Hamiltonian J . For sufficiently high spin this drastic 
simplification describes the physics at low energies very well. 

In Figure ||] we compare our predictions to the spectrum of H\ as 
obtained through numerical diagonalization. The comparison is good, 
particularly along the first excited state and near the isotropic limit. 
The reason that the Ising limit compares poorly is that the quantum 
fluctuations have the effect of regularizing the eigenstates if), whereas 
for the Ising model these states are definitely not smooth. However, 
for any q > 0, if J is made large enough, then we believe that these 
asymptotics eventually dominate. In Figure |||, we have plotted the 
spectrum of J for fifty sites and r = 0, which corresponds to the 
magnetization for which one has a minimal gap. It is clear that there 
is only a single isolated eigenvalue beneath the branch of (what would 
in infinite volumes be) continuous spectrum. Of course, to recover the 
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Inverse anisotropy : Ifl Inverse anisotropy : 1,'A 



Figure 5. Solid lines are the predicted values of spec- 
trum according to Boson gas model, circles are actual 
values of the spectrum for full XXZ as obtained by Lanc- 
zos 




0.1 02 0.3 0.4 0.5 O.S 0.7 0.8 0.9 



Figure 6. The spectrum of the Boson coupling matrix 
versus anisotropy 

spectrum of 'K^ one must take all (nonnegative) integer valued linear 
combinations of these lines, since each line corresponds to the first 
excited energy of an independent boson. In Figure [7|, we have plotted 
several multiples of the eigenvalue line, to show how many eigenvalues 
lie beneath the continuous spectrum. There is also interesting behavior 
for other values of r. In Figure |] we have plotted the spectral gap of 
J as a function of r and A -1 Finally we mention that this analysis 
can be done for any dimension, not just one. Thus one may obtain, at 
least heuristic, information about the low spectrum of quantum spin 
systems in higher dimensions by analyzing J7a, which is the first order 
quantum correction to the classical ground states for large but finite J. 
The higher-dimensional case, where the low-lying spectrum is known 



47 




Figure 7. The spectrum of the Boson coupling matrix, 
with multiples of the lowest boson energy included 

to exhibit additional structure ||, is the subject of a separate paper 
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